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Valuon  o.f  tolnl  electron  content  of  the  lonoephie 
re,  obtained  at  the  Obsorvatorio  del  Ebro  by  the 
Faraday  rotation  method,  have  been  analyzed  to 
ntudy  electron  production  rates.  Daily  values  of 
the  electron  production  rate,  integrated  with 
respect  to  height  through  the  ionosphere, were  de_ 
termined  at  sunrise  for  the  period  August  1973  - 
July  1974,  by  a method  that  makes  use  of  a two 
components  model  atmosphere,  \ Fourier  analysis 
of  the  results  indicates  an  annual  and  a semian- 
nual variation  of  the  production  rate  values.  A 
comparison  with  results  of  other  authors  seems 
to  indicate  a latitudinal  dependence  of  the  semi^ 
annual  variation;  no  similar  effect  was  found 
for  the  annual  variation. 
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1.  DITPODUC'l^IOK. 

In  a previous  report  (Aileron  and  (JaldCn,  lS;74)we  g£ 
ve  results  of  total  electron  content  of  the  ionosphe 
re  for  a period  of  six  months.  The  data  v/ere  derived 
every  2,5  minutes  from  the  signal  of  the  ceo8tation£ 
ry  satellite  IHTELSAT  II-I'3  by  the  Faraday  rotetion 
method,  /in  analysis  similar  to  the  one  described  by 
Garriot  and  Smith  (iy65),  van  applied  to  the  sunrise 
period  date  in  order  to  obtain  daily  values  of  the 
integrated  ionization  rate  of  the  ionosphere  for  an 
overhead  sun. 

In  the  present  v/ork  we  use  a more  sof  1 aticated  me- 
thod to  obtain  more  accurate  valuer,  of  the  integra- 
ted ionization  production  rate  fcr  a zenital  angle 
of  the  sun  of  9C2,  (Qqo)*  IJi  the  follr>?/ing  sections 
we  describe  the  method  of  analyoip,  in  v/hich  a two 
component  model  of  the  neutral  atmosphere  has  been 
incorporated.  The  model, whose  desciiption  is  given 
in  2,2,  la  used  in  2.3  and  2,4  to  obtain  an  expres- 
sion for  the  photoionization  rate  for  grazing  inci- 
dence. An  expression  for  the  ionization  loss  rate  is 
alao  obtained  in  2,5#  The  formulae  obtained  in  these 
paragraphs  are  used  in  2,6  to  obtain  the  theoretical 
J total  electron  content  given  by  the  model  for  the 

^ sunrise  period.  The  part  2.7  contains  the  method  to 

I obtain  the  integrated  Ionization  production  rate  by 

^ comparing  the  theoretical  and  the  experimental  data 

I of  the  total  electron  content.  V/e  give  in  section  3 

' the  results  of  the  production  rate  found  at  the  Ob- 

I oervatorlo  del  libro  by  the  application  of  the  method 

t deaorlbad*  Its  mean  value  as  well  an  its  seasonal  va- 
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riations  are  compared  with  the  reaulta  of  other  hu- 
thora  and  its  implications  are  discussed.  Finally, see 
tion  4 containo  the  concluseions  drav/n  from  the  pre- 
vious disouesion. 

2.  tIKTIIOD  OF  AlIALYSIS. 

2.1  THE  CONTINUITY  E(JUATI0N, 

The  determination  of  Q90  is  done  by  comparing  the  ex 
perimental  data  of  total  electron  content  (TEC)  dedu 
ced  through  the  Faraday  rotation  method  for  the  sun- 
rise period,  to  the  theoretical  values,  deduced  from 
a model  of  the  ionosphere , at  the  same  period. 

The  model  is  based  on  the  well-kiiown  continuity  equ£ 
tion  for  electrons: 

-^  = <l-l-V(Nvi  (1) 

were  q and  1 are  respectively  the  production  and 
loss  rate,  N the  electron  density  and  the  drift 
velocity  of  the  ionization,  produced  by  all  the  dif- 
ferent forces, 

Thio  equation  is  valid  at  any  height,  z,  of  the  io- 
nosphere, so  tiiat  it  can  be  integrated  over  all  hei- 
ghts, Before  doing  it,  we  assume  that  only  vertical 
derivatives  contribute  to  the  value  of  the  divergen- 
ce, Then,  the  divergence  term  can  be  neglected  when 
Integrated  over  all  heights.  This  in  equivalent  to 
assume  that  the  T,£,C.  Is  not  modified,  by  changes  in 
the  electronic  profile,  due  to  vertical  movements  of 
ionization.  In  the  real  ionosphere,  vertical  move- 
ments of  ionization  do  have  an  indirect  effect  on 
the  electron  content,  because  the  loss  ooeffiolent 
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dlntlnishes  with  increafiinc  hei/^hts.  Nevertheless, this 
effect  is  very  small  at  sunrise  conditions,  first  of 
all  because  the  loss  itself  is  very  small,  and  secon 
dly  beacause,  due  to  the  rapid  increase  of  the  elec- 
tron content  during  this  period,  most  of  the  ioniza- 
tion pressent  at  e given  time,  has  been  produced  lit- 
tle time  before.  If  the  movement  is  not  very  fast, the 
displacement  of  this  ionization  should  not  be  big 
enough  to  produce  appreciable  changes  on  the  loss  ra- 
te, In  fact,  at  middle  latitudes,  at  sunrise,  only  on 
E-\7  movement  of  ionization  could  affect  substantia ly 
the  total  electron  content,  because  of  the  strong  griB 
dient  of  electron  density  in  that  direction  at  that 
period.  But,  unless  there  ere  strong  electric  fields, 
that  we  do  not  expect,  this  movement  is  impeded  by 
the  geomagnetic  fiel,  whose  direction  lays  on  a plane 
almost  perpendicular  to  the  E-W  dirtction. 

The  integration  of  the  two  membres  of  the  equation(l) 
with  the  conditions  just  assvimed  gives 


dNr 

*Tf“ 


= Q-P 


(2) 


V/here  N|p  is  the  total  electron  content,  and  Q and  P, 
respectively  the  integrated  production  and  loss  rate, 

2.2  NL'UTH/lL  ATIIOSPHLIIE, 

The  two  terms  of  the  second  member  of  equation  (2)  de- 
pend very  strongly  on  1;he  properties  of  the  neutral  at 
noBXihere  that  receives  the  ionizing  radiation.  To  sol- 
ve equation  (2)  is  then  necessai'y,  to  establish  a mo- 
del of  neutral  atmosphere  from  which  expressions  for 
the  different  termes  can  be  deduced. 
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The  model  we  have  adopted  Is  an  aproxlnation  of  that 
of  Jr.cchia  (1971).  ’wo  have  found  oome  relatively  aim 
pie  formulae  that  give  o reprei-entation  of  the  neu- 
tral atmosphere,  similar  to  the  Jacchia  model  for 
the  same  boundary  conditions.  The  formulae  are  valid 
for  an  exospheric  temperature  range  at  least  between 
about  7002  K and  10002  K that  corresponds  to  the  tern 
perature  variation  during  the  considered  period. 

The  lower  limit  of  validity  of  the  model  proposed  is 
Zq  = 120Km,  The  temperature  at  this  level  is  not 
keep  constant  as  in  other  models.  On  the  contrary, 
we  make  it  dependent  of  the  exospheric  temperature, 
in  order  to  obtain  a better  fitting  to  the  experimen 
tal  data, 

A list  of  the  variables  that  are  use  on  the  formulae 
is  given  below: 

m = molecular  weight 

K = Bolzman  constant 

T,  = exospheric  temperature 

T = temperature  at  height  z 

To  = temperature  at  height  zq  (depends  on  T„  ) 
g = gravitational  acceleration 

n(X)  = numerical  density  of  the  X component  at  height 
z 

ho(x)  = numerical  density  of  the  X component  at  hel^t 
Zq  (depends  on  T»  ) 

C = constant  « 910 

b a parameter  (depends  on  T*  ) 

k » parameter  (dependa  on  T»  and  on  m) 

s *>  parameter  (dependa  on  T*  and  on  m) 

All  the  variables,  as  well  as  the  temperature  and  den- 
sity at  the  Zq  level,  can  be  obtained  from  the  exosphe 
rio  temperature,  T»  , through  the  formulae  given  in 
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pendix. 

The  exospheric  temperature,  depends  on  the  solar  radl£ 
tion  and  the  geot»gnetic  activity,  and  the  method  to 
calculate  ;-t  in  the  sajue  as  in  the  Jfecchia  racdel. 

The  temperature  distribution  with  height  is  given  by 
the  equation 

T = (3) 

were  Zq  = 12CKm,  To  and  b depend  on  T and  are  ob- 
tained thi’ough  the  equations  given  in  Appendix. 

In  fig,  la  comparison  between  the  temperature  distri- 
bution given  by  (3)  and  tlie  Jscchia  (1971)  model  is 
given  for  two  different  exospheric  temperatures.  As 
can  be  seen,  in  both  cases  the  curves  are  very  simi- 
lar. 

■Ve  ansune  the  neutral  atmosphere  to  be  composed  by 
three  elements, 0,  Later  on  v/e  shall  nee  that  the 

contribution  of  O2  to  the  total  electron  content  va- 
riation is  really  negligible  so  that  only  0 and  N2 
shall  be  taken  into  account  in  the  analysis  of  the  da 
ta,  nevertheless,  in  this  part,  the  three  elements 
are  considered. 

We  assume,  as  usual,  that  the  three  elements  are  in 
diffusive  equilibrium,  so  that  the  density  distribu- 
tion v/ith  height,  z,  of  each  component  is  given 

/'*  (4) 

fT  /Ti  dz 

n = no(V^le  4 

To  obtain  an  expression  of  the  integral  that  appears 
in  the  exponential,  v/e  aproximate  the  relation  g/T 
through  the  sum 
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g/T  rC/TW(z-zJ  (5) 

where  0 is  a constant  and  f(z  - Zq)  a function  who- 
se Integral  is: 


yilz-zoldzsln  [f|z-Zo)] 


(6) 


with 


F|z-zJ=Uk|z-zJe®‘^’^o^ 


(7) 


where, as  has  been  indicated,  k and  s are  two  parame- 
ters that  can  be  obtained  from  T„  through  the  for 
mulae  of  Appendix, 

With  these  a proximal  ions,  the  integral  of  ea_.  4 be- 
comes 

Jgn  dz  sC/(bT.)ln(T/T,e‘’‘^‘^o)),in(F(z-Zo))"'^^"' 
and  the  lensity 


{UJDQ£.)  .i22_(2-z,) 

n=nT  T e F|z-z,) 


(8) 


V/e  put 


• mC 


(9) 


I 


that  has  the  form  of  a scale  height  of  temperature 
T.  and  aoeleration  of  gravity  ecj.ual  to  0, 


Temperature  profiles  of  the  neutral  atmosphere  and  its  comparison  with  Jacchia's  (1971)  model 
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If  we  call 

and  gubatitute  P ( ?:  - zq)  by  ita  expresHion(t) , the 
equation  C be  comer. 


I 


i 

i 


I 


n =nJ 


r.-L,,. 


H. 


(z-z.) 


(Uk(z-zJe' 


.s(z-zj 


(11) 


with  T given  by  (4)» 

In  Tables  I v/e  give  the  values  of  log  n for  the  three 
oonponents,  obtained  by  eq,  11  for  different  exoaph£ 
Tin  temperatures.  /'Iso  the  corresponding  values  of 
the  Jacchia  (1971)  model  and  the  relative  difference 
between  both  models  are  shown,  i.s  can  be  seen  the  dif 
ferences  between  the  models  are  not  greater  than  2^5 
till  very  high  altitudes,  and  the:'  are  greater  than 
lOfJ  only  in  regions  were  the  density  is  less  than 
20crQ,  The  contribution  of  this  region  to  the  total 
ionisation  of  the  ionosphere  is  so  small  that,  even 
if  the  errors  in  the  determination  of  the  density  of 
the  componens  are  very  high,  the  results  are  not  af- 
fected. 


2.3  PHOTOIOhIZATIuN, 

AS  it  ia  well  known,  most  of  the  free  electrons  in 
the  middle  latitude  ionosphere,  are  produced  by  pho- 
toioni7.ation. 

In  a mixture  of  gases,  the  photoionisation  rate  of 
each  component,  qj  , produced  by  a monochromatic  ra 
diation  of  flxnc  ^ is  given  by 


T 
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q =A  a.n.e 
J • ‘J  I 


-X.a  /n 

(ai4 


ds 


(12) 


where  a—  and  ^ are  respec  tively  the  ionization  and 
abaortion  croas  sections  of  t'le  j component  and  nj 
the  numerical  density.  The  sun^l  ia  extended  to  all 
the  components  and  the  integral  is  taken  along  the 
ray  path. 

Since  the  ionizing  radiation  of  the  atmosphere  is  not 
monochrojnatic,  eq,  12  cannot  be  applied  directy,  be- 
oaune  the  cross  sections  depend  on  the  wave  length, To 
use  that  equation,  v;e  have  found  an  equivalent  cross 
section  for  the  total  ionizixig  flux.  The  ionizing; 
spectrum  has  been  divided  into  intervals  of  narrow 
bandwidth,  so  that  the  cross  sections  can  be  conside- 
red constant  inside  them.  The  cross  section  data  co- 
rresponding to  each  interval,  a|^  , as  well  as  the  so 
lar  radiation  flux  for  the  same  frequecles,  , ha- 
ve been  taken  from  Kookarts  (1973)  and  a weighted 
mean,  e , of  the  cross  sections  has  been  taken  for 
eacli  element  so  tliat 


where  is  the  total  ionizing  solar  flux.  Then 
equation  12  can  be  used  taking  s as  crons  section 
and  as  the  ionizing  flux.  The  values  of  the  absor 
tion  cross  sections  found  for  the  three  elements 
are  t 


aQ=6.17 10'’®crTv^ 


=11.510‘’®cmr^ 

N2 


or  =13.3910"’®cm"2 
02 


’’K 


I 

1 

I 

i 

I 
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The  values  for  0 and  ’^2  higher  than  those  of 
Smith  (1968) 

(To  = 5 .7  10"’®  cm'^  ; a , -7.4  10"’® cm"^ 
u ' N2 

and  lov/er  than  those  of  Yeh  et  el  (l^fio) 

OQ=7.3210r  cm  ; o^^r14.110  cm 


2.4  GlhlZINO  INUlDi;iiC9. 

2.4.1  'bsortion  of  radiation. 

The  optical  deptii,  ■ts<Tq/^nds  , in  eq.  12  gives 

the  ray  abaortion  along^Vne  path. 

Since  we  have  tr'  obtain  t’le  ionization  at  the  sun- 
rise period, grazing  inci.dence  of  the  radiation  has 
to  be  considered.  V/e  shall  assume  tiiat  the  ray  path 
is  a straight  line. 

In  fig, 2,  0 represents  the  centre  of  the  earth, 0-S 
the  line  joinning  the  centre  of  the  earth  with  the 
sun,  IITII  the  surface  of  the  earth  and  It  the  earth  ra 
dius,  V/e  shall  calculate  the  absortion  of  the  ray 
till  the  point  ?,  To  do  it,  we  shall  find  an  exjtres 
Sion  of  tlie  density  in  a point  P'  of  the  ray  path, 
that  can  be  integrated  along  the  ray  path. 

From  fig,  2 

prIzVR  )sin  A R)sin  "X 

if  we  call 

h = 4ili  : h..-.a±24  (13) 

H.  ' * H, 

we  obtain 


I 


Grazing  incidence  geometry 


12  - 


liSr 'I' -M 

I’o  find  the  density  in  ?'  it  ia  enoutsh  to  substitute 
- !?o  in  tho  different  termo  of  the  expression  of 
n (ofr.  11).  For  'i’j^  we  obtain 


,j\  _ Sin> 

T. 


-b|z-z,l 

e 


-r 


The  term  to  be  subtracted  in  this  equation  is  always 
smaller  than  1 (and  much  smaller  when  the  altitud 
is  higher)  so  that  only  the  first  terms  have  to  be 
retained  in  the  development  of  the  power  and  the 
equation  becoiMS 


r'=T;'(i.ri^ 


e 


sinX 

sTnT 


-1) 


-t)(Z-2,) 

e ) 


The  substitution  of  z'-'/.q  in  the  other  terms  of  eq, 
11  Is  straightf on/ardy  and  doing 


-j-*b  = B ; -s=S  ; -i-  ♦b-ssG 

n«  H.  H, 


= T. 


It  gives 


13 


n,.=nT’^f’'e  e H.*j 


4^  -Blz-z.)  -BA  -S(2-z.)  -SA 

e ♦k(z-z,)e  e ♦ 


(14) 


-G(z-z,)  -GA  -SIz-z,)  -SA  -&|z-z,)  -GA 

R 


« **<6  Ae  *10  ke  Ae 


The  ds  that  appears  in  T is  (cf,  fip,  2) 


ds  = <R*z)sinX 

sin  \ sin2x 


or,  acGordinc;  to  (13) 

ds=  H.hsinXcosec^A  dx 

v/ith  the  limits  of  the  integral  of  x O^A^X 


To  calculate  the  integral  we  ascu,Tie  that  is  in 

dependent  of  A and  its  value  is  equal  to  the  one 
that  corresponds  to  A = X • Since  R $>  z in  the  re 
gion  where  most  of  the  ionization  is  accumulated, the 
Slim  R ^ z is  taken  as  constani;  and  equal  to  6,700 
Km,  With  these  conditions,  we  find  the  following  ex- 
pression for  the  optical  depth 


-IzIm.  -BIz-z.)  -Ste-z.) 

'r=K.K,e  H,  ♦Kje  ♦K3(z-z,)e 


(15) 


-G(z-z,)  -S(z-z,)  -G(z-z,) 

♦K^(z-z,)e  ^*^5*  ^*^6* 
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where 


rr-T 


KQ  = H,hsinx 

Kj  =/  e cosec^X  dX 


cosec 


2x  dX 


K j =k/  e cosec^  A dA 


1 

K^rT^k/ 


K5=k/ 

'a 


-GA  ^ 

e cosec^X  dX 


SA 


Ae  cosec^XdX 


i' 


r*  A 

Ke  sTr  k/  Ae'  cosec^x  dX 


All  these  Integrals  arc  independent  of  z and  can  be 
obtained  by  a nximorical  integration.  The  Integrals 
to  K4.  have  the  same  forru  and  only  differ  in  the 
constant  coefficient  that  multiply  the  factor 

The  products  Kq  Ki  ( i=l,i)  have  the  form  of  the  Cha£ 
man  function  (efr.  Chapman  19^1).  In  fact,  Kg  is 
the  Chapman  function  for  a scale  height  H.  , This 
means  that  integraln  Kp  to  give  the  modifica- 

tion of  present  model  as  compare  ’.vith  0 Clxapm  .n  1l- 
yer. 


T 


r 
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2.4.2 


groduot.i.on  rate. 

The  expresHion  for  the  production  rate  is  easily  ob- 
tained by  oubstitutinj!^  the  optical  depth  of  eq,  15, 
into  eq.  12.  For  the  production  rate  of  the  j com- 
ponent it  given 


-b|z-z,) 

-I,K.je  ”-lK,j*K2je  ♦(K3|(z-z,)* 


(17) 


stz-j;) 


♦Kgjle  ♦IK4||z-z,hKgj)e 


-(b-s,Mz-z,) 


were  nj  is  a function  of  height  and  is  given  by(ll). 

To  compare  the  ionization  production  rate  of  the 
three  coraponontn, 0,  1^2  and  Opr  the  distribution  cur- 
ves of  qj  with  height  have  been  obtained  for  seve- 
ral exospheric  temperatures  asauriing  an  ionization  ef 
ficioncy  equal  to  1.  The  results  are  shown  in  fig. 3. 
For  the  solar  flux,  the  value  0=  5.24 10*®  em'^s"' 
has  been  taken,  that  corresponds  to  the  sun  of  the 
fluxes  for  the  different  wave  lengths  between  910  and 
00  n given  be  Kokarts  (1973). 

As  can  be  seen, the  contribution  of  the  Op  production 
rate  to  the  total  ionization  is  very  small, so  that  it 
can  be  neglected  for  the  calculation  of  the  total 
electron  production. 

For  such  a calculation  we  liave  assumed  that  the  ioni- 
zation efficiency  of  the  atomic  oxigen  is  equal  to  un^ 
ty  and  tliat  of  Np  equal  to  zero,  that  is 


t 
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The  firtit  hipotheeis  agrees  with  the  experimental 
suite  (cf,  hhitten  ami  Poppoff  (1971).  The  second 
one  ie  an  aproximatior.  based,  first  of  all, on  the 
fact  that  only  scuill  nuiiber  of  ^2^  ions  Jiave  been 
found  in  the  thermosphere  (cfr,  for  instance  Narcisi 
(197!;)).  This  indicates  that  N2^  ions  disappear  so 
quickly  that  its  contribution  to  the  total  ioniza- 
tion is  negligible.  On  the  other  hand  the  reactions 
for  the  N2  ions  diseppearance,  do  not  produce  other 
long  life  ions  that  could  contribute  to  the  ioniza- 
tion. 

With  these  assumptions,  the  production  rate  of  elec- 
trons is  practically  the  seme  as  the  ionization  pro- 
duction rate  of  the  atomic  oxigen.  Its  expression  is 
the  sane  of  eq.  17  with  the  j component  being  the 
atomic  oxigen  and  the  sum  of  the  exponential  being 
extended  to  0 and  H2, 

In  order  to  compare  with  the  experimental  data  we 
need  the  total  electron  content.  The  production  rate, 
q , has  then,  to  be  integrj«ted,  to  ob-tain  the  inte- 
grated production  rate  Q.  This  integration  yields 


W /r**  ^ 


(10) 

dz 
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vfhere  v/e  put  Qx  to  indicate  tho  dppend'?nce  of  the 
zenitel  angle  of  the  sun  through  end  v/here  all 
dependence  of  z hoc  been  explicited. 

In  fig.  /*,  the  curvec  Q ( X ) between  852^X<100e 
are  cho?;n  for  tv/o  exocpheric  tempers tureo,  Theee 
curves  are  obtained  from  eg  18  considering  the  up- 
per limit  of  integration  x..=  500Kin  since,  as  it 
can  be  seen  in  the  curves  of  fig.  5»  the  production 
rate  above  this  region  is  almost  zero. 

The  two  curves  of  fig.  4 are  almost  parallel.  This 
means  that  both  curves  can  be  reduced  to  a common 
one  if  they  are  normalized  by  a convenient  factor. 
Smith  (196b),  with  a model  of  an  isothermal  neutral 
atmosphere,  with  two  components  whose  production 
profiles  are  derived  from  the  Chapman  theory,  fovuid 
that,  normalizing  the  Q curves  with  respect  to 
its  value  for  *)(  = 902  , the  shape  of  the  curves 
Q / O90  almost  independent  of  the  concentration 
relation  of  both  component n,  V/e  have  applied  a sim^ 
lar  normalization  to  the  present  model  and  have 
found  that  the  curves  Q / Qqc  sre  almost  independent 
of  the  exospheric  temperatvire.  This  cm  be  seen  in 
fig,  5 where  the  curves  corresponding  to  exospheric 
temperatures  T®  = 700‘’K  and  Tu?  =10002  K are  shown. 
This  property  of  the  quasi- independence  of  Q / Q^q 
from  the  exospheric  temperature,  shall  be  utilized 
later  on  to  deduce  the  value  of  Qgo* 

.5  lOHIZATIOH  LOSS  RATE, 

As  it  la  known  moot  of  the  0'*'  ions 
pear  throvigh  the  process  of  reactions 
0%N2  — N0**N  a) 

N0%t-—  N*0  b) 


deaap- 

(19) 


fig.  5.  - Integrated  ionization  production  rate,  normalized 
with  respect  to  Q^q.  as  a function  of  the  zenital 
angle  of  the  sun 
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llie ion  0'*'  ie  dominant  in  the  recion  between  about 
200  and  500Kti,  v/here  most  of  the  ioni?:ation  is  pre- 
sent, It  eeems  then  reeeonably  to  conclude  that  reac 
tion  19-b  is  much  more  rapid  than  reaction  19-3 » eo 
that,  as  a first  aproxiniation,  we  can  assume  that 
the  ionization  loss  rate  is  eauol  to  the  rate  of  oc- 
currence of  reaction  19-b» 

According  to  this,  the  ionization  loss  rate  is 

l=CQn[o+]n[N2]  (20) 

where  Cg  is  the  loss  coefficient  of  reaction  19-a 
and  n [ X ] indicates  concentration  of  the  X ele- 
ment. 

As  we  have  said,  cq,  20  is  only  an  approximation  but, 
since  during,  the  sunrise  period  the  ionization  loss 
is  very  small,  an  error  in  its  determination  cannot 
be  very  important. 

Because  of  the  neutrality  of  the  ionosphere,  we  can 
also  approximate  n[0'*^]  n[e~  ] and  eq,  20  beco 
mes 

l = nn[e-]  (21) 

where 

0 = Car[N2] 

For  the  variation  of  f)  with  height,  we  have  used  the 
simple  model 

- z-z. 

flrO.e  H,n2  (22) 

whore  scale  height  of  112  correspon- 

ding to  the  temperature  To*  , The  use  of  the  varia- 
tion of  the  H2  density  given  by  the  more  sofistl- 
oated  model  of  this  report  to  calculate  n , would 
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ooraplicate  the  calculation  without  apreciably  ircre^ 
asing  its  accuracy. 

The  value  of  the  coefficient  ftv  , that  appear  in  eq. 
22,  has  been  obteine^  from  the  value  of  at  a height 
of  30oKdi,  given  by  Smith  (1968),  Oitt  ^ 

2,6  TOTAL  LLliCThOK  COIITINT  VjlRWI'IOU, 

With  the  simplification  related  to  tlie  movement  term, 
and  taking  Into  account  eq,  21  am'  22,  the  continuity 
equation  (1)  becomes 

N(z.t)  (23) 

whore  the  dependence  from  the  z and  t variable  has 
been  shovm#  The  expression  for  q (z,t)  is  given  in 
eq,  17  and  dependens  on  t through  the  zenital  angle 
of  the  sun. 

To  compare  with  $he  experimental  total  electron  con- 
tent data,  we  have  to  solve  the  continuity  equation 
(23)  and  calculate  afterwards,  the  integral 
for  the  different  values  of  t at  which  data  have 
been  recorded. 

To  solve  eq,  23,  we  fix  as  the  initial  moment,  to,the 
time  when  X =1002,  since  the  photoionization  is  prac 
tically  zero  for  greater  values  of  x and  the  experi- 
mental total  electron  content  usually  do  not  increase 
before  this  time. 

Using  a method  similar  to  the  one  utilized  by  Smith 
(1968)  the  total  electron  content  is  divided  into  two 
parts,  the  nighttime  ionization,  or  ionization  present 
at  the  initial  moment  t=  to,  and  the  loAlzatlon  pro- 
duced from  this  moment  onwards.  Both  parts  can  be 
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treated  separately  becauee  of  the  linearity  oi  eq, 

23. 

The  second  part  of  the  electron  content  can  be  found 
by  nuaierloally  solving  eq,  23  with  the  initial  oond^ 
tion  N(z,to)  = 0,  and  integrating  afterv/ards  N(z,t) 
with  respect  to  z.  For  the  nighttime  ionization,  we 
asBune  that  its  evolution  during  the  time  between 
1002>  X>  872  (at  which  the  cplculations  are  exten- 
ded) is  the  8amt>  as  during  the  hour  before  the  ini- 
tial moment  t=  to  {ot  \=  1002 ),  The  limit  of  y,  =87® 
is  Imposed  because  the  validity  of  tlie  aproxlmatlon 
(21)  la  questionable  afte].’  that  time.  In  the  fig. 6 
an  exeraple  of  the  variation  of  Ti;o  from  an  hour  be- 
fore “X  *1002  is  shown.  As  it  can  be  seen,  the  vari£ 
tion  till  y.  =1002  can  be  fitted,  by  a straight  line 
of  negative  slope.  If  no  new  ionization  where  produ 
oed,  the  electron  content,  Kp  , at  any  time,  t , of 
the  Interval  100 X 872  should  be 

Np  = N^^ 

where  Njq  is  tha  electron  content  at  the  initial 
moment,  tot  and  -b  ic  tlie  slope  of  the  fitted  li- 
ne. Hp  is,  then, the  amount  of  ionization  present 
at  time  t , due  to  the  nighttime  ionization;  adding 
it  to  the  ionization  produced  from  the  Initial  mo- 
ment, tQ,  till  time  t , we  obtain  the  total  elec- 
tron content  at  this  time. 

The  slope  -b  is  obtained  by  fitting  a straight  line 
to  the  experimental  TEC  data.  can  be  obtained 

rectly  from  the  records  or  by  the  procedure  that  we 
give  below. 
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fig.  6.  - Total  electron  content  (Nj)  and  its  values  corrected  for 
the  effect  of  the  nighttime  ionization 
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If  we  aubntlfcute  q (z.t)  by  q(z.t)/  Q90,  equation 
2!>  becomes 


d NrIz.O 


q(z.t) 

Q90 


-ne 


z-z. 

^•N2  Nplz.t) 


(24) 


with 


N^lz.t) 


N(z.t) 

Q90 


Solvlntj  this  equation  v/lth  the  initial  condition 
Nr  (Kjto)  = 0,  instead  of  Nj  (t)  we  obtain 

N,r--Nt,«)/Q9o 

where  Njj  indicates  the  total  electron  content  pr£ 
duced  by  photoionization  from  the  initial  moment  t© 
till  time  t.  Then,  the  total  electron  content  is: 

Nj(t)  = N^^.b(t-t.)4Q9oNjj^(t) 

and  putting 

NTc<t)rNy(t)4b(t-t.)  (25) 


results 

NTc(t)=NvQ90Njj^(t)  (26) 

that  gives  Qgo  as  the  slope  of  the  linear  correla- 
tion between  the  values  of  and  UTr* 

The  values  of  are  obtained  by  oorreoting  the 

experimental  data  of  TEC  for  the  nocturnal  lonlza- 
tion«  as  it  is  nho^Tn  in  fig.  6.  The  are  the  the 

oretloal  values  of  TEC  normalized  with  respect  to 
QgOf  obtained  by  solving  eq.  24'. 

The  correlation  of  these  series  of  values  gives  also 
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the  TEC  at  the  initial  moment  NipQ, 

The  use  of  eq,  24  instead  of  eq,  23  han  the  advan- 
tage that  the  values  of  Njp  depend  on  the  0/Qgo 
relation  that*  as  we  have  indicated, is  almost  lnd£ 
pendant  of  the  exospheric  temperature,  T«  , from 
which  the  parameters  of  the  equation  are  derived. 

On  the  other  hand,  in  eq,  24,  q values  are  subs- 
tituted by  q/Qgo  and,  while  q directy  depends 
on  the  solar  flux,  q/Qgo  is  independent  of  this  pa^ 
raraeter.  Also  the  dependence  of  q on  the  product 
of  the  oxigen  cross  section  and  the  oxigen  density 
at  the  Zq  level  ( tfo*  ‘ ) is  stronger  than  de- 

pendence of  q/Qgo  on  the  same  product,  There-fore, 
any  error  In  the  determination  of  these  parameters 
has  much  less  influence  in  the  value  of  Qqq  If  is 
calculated  by  eq,  24  that  if  it  is  done  by  eq.23, 

3.  RESULTS. 

3,1  VALUES  AT  0B3ERVAT0RI0  Dili  EBRO, 

The  method  Just  described  has  been  applied  to  the 
TEC  data  obtained  in  the  Observe torio  del  Ebro  du- 
ring the  period  Au/^iust  1973-July  1974,  to  deduce 
dally  values  of  Qgo.  As  we  have  seen,  this  requires 
the  determination  of  daily  values  of  during  the 
interval  lOOB  ^ 87^  which,  in  turn,  requires 

the  determination  of  the  exospheric  temperature 
for  the  same  interval.  As  we  have  already  said,  T« 
has  been  oalculated  by  the  method  given  by  Jaochia 
(1971),  In  this  method,  it  is  assumed  that  there  is 
no  semiannual  variation  of  the  exospheric  tempera tu 
re.  On  the  contrary,  in  his  model  of  1970,  Jaochia 
Incorporates  a aemlannual  variation  of  Tm  to  ao- 
count  for  the  semiannual  variation  of  density  dedu- 
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ced  from  the  satellite  drag  data.  The  semianrual  va- 
riation of  T«,  , calculated  according  to  the  Jacchia 
(1970)  model,  gives  a maximum  oscilation  of  t 409 
for  the  ueriod  August  1973-  JUI3'  1974.  The  incorpora 
tion  of  this  variation  to  T«»  would  not  change  Bub£ 
tantlal]y  the  rcaulta  obtained  in  this  report,  since, 
as  has  been  shovm,  the  influence  of  Ta>  on  the  Q^q 
calcf.lation  is  very  small  in  the  method  used.  Per  the 
same  reason,  we  have  neglected  the  small  variation  of 
Ta>  during  the  interval  lOO^  X 872  and  have  taken 
it  ne  constant  equal  to  its  value  at  X = 90e. 

In  fig. 7 some  examples  of  the  representation  of 
the  Njq  and  values  are  given.  The  slope  of  these 
linec  ia  the  value  of  Qgc*  In  some  cases  it  seems  that 
Nipj^  only  startD  increasing  several  minuter  after  =1008 
In  these  cases  two  linear  correlations  have  been  obtai 
ned,  the  first  one  including  all  the  points  while  in 
the  second  one  the  points  before  the  first  four  Increa^ 
sing  values,  have  been  eliminated.  In  general  both  co- 
rrelations are  very  similar  and  only  in  few  cases  where 
the  first  one  gives  an  absurd  result  ( !q(y<0)  the  dif- 
ference has  been  significant.  Therefore,  the  second  co 
rr»*latton  has  been  taken  as  the  valid  one  in  determi- 
ning Qgo. 

The  values  of  (jgo  obtained  for  the  i/holc  period  August 
1973-July  197^  are  shown  in  fig.  8 wliere  also  the  run 
ning  mean  of  27  days  of  these  values  is  given.  The  sm^ 
othin  , of  th?  QgQ  values  has  been  used  to  eliminate 
the  possible  influence  of  the  rotation  period  of  the  sun. 
Although  there  is  a certain  dispersal  of  the  points,  ao 
me  sort  of  semiannual  variation  is  clearly  appreciated 
tliat  is  bettor  seen  in  the  nmning  mean  curve.  In  this 
curve,  one  maximum  oooura  about  the  20  of  October  and 
another  one  at  the  end  of  Pebroary  or  beginning  of 
Uaroh.  This  aooond  one  seems  to  be  a little  lower.  The 


es  and  its  Z7  days  running  mea 
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minimum  at  the  end  of  June  or  the  beginning:  of  July 
is  also  lower  that  the  one  at  about  12  of  January, 

The  minimum  of  oumnier  is  leoo  sharp  than  the  one  in 
winter  perhaps  because  of  a bigger  dispersal  of  the 
values  since  the  month  of  March.  The  dispersal  may  be 
caused  by  reflections  of  the  incoming  wave  in  an  obs- 
tacle placed  towards  the  west  of  the  receiving  ante- 
nna, The  satellite  had  a shift  towards  the  west  and 
we  have  to  move  the  antenna  in  that  direction,  so  that 
the  values  of  the  last  months  can  be  a little  affected 
by  reflections  in  the  mentioned  obstacle. 

If  the  difference  between  the  values  of  suriraer  ana 
winter  is  significant,  it  could  indicate  the  existence 
of  an  annual  variation,  besides  the  semiannual  one, al- 
though smaller  than  that, 

iOi  harmonic  analysis  of  all  the  points  gives  a mean 
value  of  Q90  = 1,010*^  effl^  a first  harmonic  of  am- 
plitude 15f»  of  the  mean  value  and  a second  harmonic  of 
20^  amplitude.  If  the  two  harmonics  are  expressed  in 
the  form  R cos  (k<«*t  the  phases  are  <^<=29«69 

and  <(^1-  185.3®  beginning  on  January  Ist,  that  corres- 
ponds to  an  ann\xal  variation  with  maximum  on  December 
lot  and  minimuM  on  July  Ist,  and  a aeraiannual  varia- 
tion with  maxima  on  28  September  and  29-30  March, and 
minima  on  28  December  and  2^'  June, 

In  fig,  9 the  running  mean  and  the  curve  resulting  of 
the  two  harmonies  are  represented  for  oomparlaon, 

3,2  COLIPARISOH  WITH  OTHER  AUTHORS. 

3.2,1  Mean  value 

Smith  (1968)  analyisea  the  Faraday  rotation  data  obtal 
ned  in  Hawaii  (21'-’H)  with  a model  of  Isotermal  neutral 
atmosphere  of  tv/o  components  (N2  and  0),  based  on  data 
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of  the  Clra  (1965)  model,  and  whose  production  prof^ 
le  he  deduces  from  the  ChapTnan  theory.  He  findo  a 
mean  value  of  Qgo  = 1»7  10^  c5^  for  the  period 
Sept,  1964-Jiuguot  1966  that  includes  a minimum  of  so 
lar  cycle.  With  the  same  method,  Spurlin^^  (1972)find8 
daily  values  of  Qgo  in  Bribie  Island  (2795  153fiE) 
for  the  period  0«t69-Sept,197O  of  high  solar  activity. 
The  days  riuining  mean  of  these  values  remains  bet- 
ween 2,10  and  3»3*10  cm^  5^,  The  differences  of  the- 
se  values  with  the  value  1,0  10^  cm'^  s obtained  at 
Ebro,  oan  be  due  to  a latitudinal  variation  of  Qgo 
that  diminishes  with  increasing  latitudes  and  to  a 
variation  v/ith  the  solar  cycle, 

Tyagi  and  Ultra  (1970),  v/ith  data  of  the  Brl-C  satell^ 
te,  find  for  Delhi  (20,621i)  sununer  Qgo  values  from 
0,9310^  for  very  1o\t  solar  activity  to  2,8  I0^om“  s^ 
for  hl^  solar  activity.  These  authors  assume  an  Iso- 
therraul  atmosphere  of  only  one  component  and  do  not 
consider  the  nigh  time  ionization,  so  that  the  values 
of  Qgo  should  probably  be  a little  higher.  Assuming 
a 30yJ  Ions  due  to  nightiue  ionization,  that  ia  about 
the  one  calculated  in  this  report,  the  minimum  value 

G PI 

of  Qgo  should  be  l,3«10'cra  s slightly  higher  than 
the  value  of  obtained  at  Ebro;  these  results  are  in 
agreement  with  the  latitudinal  variation. 

Other  authors,  employing  in  general,  an  Isothermal 
neutral  atmosphere  of  only  one  component,  find  the  in 
tegrated  production  rate  for  an  overhead  sun  (Qo),  in 

order  to  compare  their  results  with  the  ones  at  Ebro, 
we  oaloulated  the  value  of  Qq  that  corresponds  to  our 
mean  value  of  Qgg  for  an  exospheric  temperature 
“ 7709K,  that  la  the  mean  temperature  during  the 
period  covered  by  our  data.  The  obtained  value  is 
Q(^  0,97  lO^^offi^  o^,  slightly  higher  that  the  one 
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lO  •X 

found  by  Titheridge  (19G6)  0,  1,10  s and  lower 

than  the  valueo  of  Garriot  and  ymith 

1.4  lO^^^oiu^  8^,  Tcj-'icr  (106F)  1,0  lO^^cSi^  c^for  thir 

teon  sumnor  days  and  2*4  li^cn^  c^for  eJ^ht  winter 

days,  laajid  and  hhuri'Tnla  (197G)  1,7  10,  1,0  10^^ 

and  1,45  10^^  cm^  respectively  for  winter,  summer 

and  equinoxes  during  a minirauLi  of  solar  activity  and 

2.7  10^°,  1.0?  lO^^and  2./1  lo^^’ciTi^  for  the  same 

seasons  and  high  solar  activity,  and  Koster  (1976) 

between  about  2,6  10  cin'^  a for  lev?  solar  activity, 

and  6 lO^^crii^  s^for  a maximum  of  the  solar  cycle. 

The  lowest  value  of  Qq  mentioned  above,  is  that  of 

Titheridge  (1966),  and  oorreeponds  to  a station  of 
the  southern  hemisphere,  nuckland  (subion,  342 o)  du 
ring  a period  of  lo-.  solar  activity,  '^ll  other  values 
belong  to  the  northern  hemisphere  and,  except  for  the 
data  of  Taylor  (1965) » correopend  to  latitudes  well 
belovr  that  of  Ebro;  Koster  (1976)  at  Legon  (526SN), 
Garriot  and  Emitli  (1965)  at  Hawaii  ( subion, 20CN)  and 
kiajid  and  Bhuriwala  (1976)  at  Karachi  (252K).  The  va- 
lues given  by  Taylor  (1966), although  correspond  to  a 
latitude  of  48211  greeter  than  the  one  at  Ebro,  they 
were  obtainei;/ during  a period  of  high  solar  activity. 
This  explains  that  his  values  are  bigger  than  the 
Ebro's,  of  lower  latitud,  but  obtained  in  a period  of 
lov;  solar  activity, 

Rao  (1967)  obtains  the  production  rate,  Qq,  for  an 
overhead  sun,  in  the  Eiaximum  of  the  F layer,  from  TEC 
data  derived  from  the  the  signal  of  the  BE-C  satelli- 
te, The  satellite  signals  v/ere  recorded  in  Urbana(40CN) 
dtxring  the  period  Auguat-Se pte ii.be r 1965  of  low  solar 
activity.  He  finds  = 0,638  10^  s^.  Assuming  a 

Ohiipman  layer  with  only  one  component,  the  production 
rate  In  the  maxisiiun  of  the  layer,  , ean  be  related 
with  the  Integrated  production  rate  for  an  overhead 
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8VUI,  Qy,  end  the  scale  heigh,  K , through  the  equa- 
tion gQ  = o H qQ  , were  e = ?.718,  Arcuniing  a value 
of  acale  heigh  H = 50Kt!i,  correnponding  to  the  oxigen, 
we  obtain  Qq  = 0,8?  10^^^  cEi.^  , wery  near  to  the  va- 
Ixie  of  Ebro  (of  aiinilar  latitude).  The  similarity  is 
higher  if  we  consider  that  in  the  Qq  calculation  of  Ur 
bana  the  nightime  ionization  luis  not  been  taken  into 
account, 

3*2,2  Annual  and  Gemicnnual  variation. 

>unong  all  the  authors  Just  mentioned,  only  Smith  (1968) 
and  Spurlin^i  (1972)  employ  a model  of  tv.o  components  to 
deduce  the  ionization  production  rate, and  neither  of 
them  find  significant  periodic  variations  of  Q , With 
respect  to  the  rest  of  the  authors,  we  have  already  men 
tinned  that  Taylor  (1969)  finds  higher  values  in  winter 
than  in  summer,  but  the  number  of  days  analized  is  so 
small  that  it  is  difficult  to  draw  any  final  conclu- 
sion, Also  Majid  and  Ehuriv/ala  (1976)  find  higher  va- 
lues in  winter  tan  in  stunuier,  independently  of  the  so- 
lar activity.  They  find  that  at  any  level  of  solar  ec- 
tiviy  the  Q values  in  equinoxes  are  lower  than  in 
winter  and  higher  than  in  sumi-'ier,  what  indicates  an  an- 
nual variation  with  maximui'  in  winter  and  minimun  In 
suinner.  At  the  came  resxilt  arrive  Risbeth  and  Setty 
(1961)  considering  the  electrcnio  density  variation 
during  sunrise  at  fix  heights,  obtained  with  a lono- 
sonde  of  vertical  incidence. 

Results  more  similar  to  ours  are  obtained  by  Titherid- 
ee  (1974),  This  author  finds,  for  a period  of  increa- 
sing solar  Hctivity,  a semiannual  variation  of  Qp  with 
maxima  in  Llarch-Aprll  and  lieptenber-f'ctober,  superim- 
posed to  an  annual  variation  v/ith  nuiximum  in  January- 
February,  for  the  stations  of  Hawaii  and  Stanford  (sub 
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ion  348N),  For  the  southern  station  of  Auckland,  he 
finds  only  the  semiannual  variation.  Although  this  au 
thor  finds  values  of  Qq  and  not  of  QcQt  the  varia- 
tion of  both  parameters  should  be  similar, 

Koeter  (1976)  performs  an  harmonic  analysis  of  the  d£ 
ta  given  by  Titheridge  and  those  obtained  by  himself 
in  the  equatorial  station  of  Legon,  Ke  finds  a first 
harmonic  (annual  variation)  of  amplitude  betrveen  12 
and  22?'-  of  the  mean  value  in  the  northern  stations 
and  only  of  3;^  in  Auckland.  The  nmplitud ; of  15:'  found 
for  kbro,  agrees  wit'n  the  values  of  t’le  northern  sta- 
tions. For  the  second  harmonic  (semiannual  variation) 
he  finds  values  betv/cen  8 and  10  for  the  four  sta- 
tions, about  half  of  the  20;'  atiplitude  of  Fbro,  With 
respect  to  the  phases  he  finds  for  the  first  harmonic 
2338  for  Auckland  and  betv/een  -142  and  19^  for  the 
others.  The  Jibro  result,  29.62,  is  still  comparable 
with  the  northern  stations.  The  phases  of  the  second 
harmonic  for  the  northern  stations  are  very  similar, 

161  and  162®,  reacimig  1722  in  .lucklan,  that  is  the 
nearef3t  value  to  the  I852  oi'  Ibro.  ^‘Iso  'uckland  is 
the  only  station,  among  the  four  analyzed  by  Koster, 
whose  semiannual  variation  is  bigger  than  the  annuel 
one,  the  same  that  happens  in  Ibro, 

Albcrca  and  Gald5n  (1974),  with  a model  of  only  one 
component,  obtained  daily  values  of  Qq,  by  fitting  a 
second  degree  polynomial  to  the  T.E.  data  at  sunr^ 

Be  at  the  Observatorio  del  hhvo.  The  variation  of  Qq 
for  the  first  six  months  of  197^  was  different  from  the 
one  found  in  the  present  report,  A more  recent  study 
of  the  method  seems  to  indicate  that  the  second  de- 
gree polynomial  apmximation,  did  not  express  the  va- 
riation of  the  TkO  at  sunrise  with  enough  accuracy  to 
obtain  the  value  of  Qq  through  the  procedure  there  in- 
dicated. In  fact,  the  coefficient  of  the  second  degree  j 

i 

I 
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term  that  was  used  in  the  determination  of  was  ve 
ry  sensitive  to  the  interval  in  v/hich  the  fitting;  was 
perforned.  In  some  cases  the  Inolusion  (or  exclusion) 
of  only  one  value  of  T'JC  was  enought  to  produce  a 
change  of  the  coefficient.  On  tlie  other  hand,  the  first 
value  taken  for  the  fitting,  was  thn  nearest  to  the  mo- 
ment at  v/hich  ')^  = 1005,  Bat  many  times, at  this  point, 
the  TIIC  is  still  decreasing,  so  that  the  coefficient 
of  the  first  degree  term  of  the  polynomial  is  negati- 
ve, As  a consequence  the  coefficient  of  the  second  d£ 
gree  term  has  to  be  bigger  to  be  able  to  fit  the  xn 
creasing  part  pf  the  curve.  This  negative  coefficient 
has  a seasonal  variation  and  its  effect  is  different 
on  the  different  seasons. 

As  an  indication  of  this  seasonal  variation,  the  slo- 
pes of  the  streight  lines  fitted  to  the  values  of  TEC 
during  the  hour  before  to  are  given  in  fig, 

10,  AS  can  be  seen,  in  April  and  Way,  the  slope  is  mo 
re  negative  than  l)i  winter,  a result  that  can  Justify 
the  form  of  the  curve  found  by  the  mentioned  method. 

If,  the  values  of  TKC  are  corrected  for  the  nighttime 
ionization  (through  a similar  procedure  to  the  one 
adopted  in  this  report) before  the  fitting  of  the  pol^ 
nomial  then,  different  valutjs  of  Qq  arc  obtained.  The 
results  obtained  by  this  metho’,  are  a little  more  sJL 
milar  to  those  of  the  present  report,  Thic  seems  to 
indicate  that,  although  the  electron  loss  at  sunrise 
is  very  small,  it  affects  the  shape  of  the  electron 
content  variation  during  this  period,  so  that  not  al- 
ways can  be  accurately  represented  by  e second  degree 
polynomial. 

Another  cause  for  the  discrepancy  of  the  results  found 
by  both  methods,  may  be  the  already  mentioned  when 


explain  the  dispersion  of  the  values  of  the  last 
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analyzed  monthu.  That  Is,  the  distoi'tion  of  tha  ele£ 
tron  content  ovirve,  duo  to  ref  lections  of  the  incom- 
raln,"  eipnal  by  an  obstacle  not  far  from  the  receiving 
antenna,  and  that  was  more  effective  during  the  last 
months.  The  cui’ves  of  fig,  11  seem  to  confirm  this  ex 
planation.  They  are  the  running  mean  curves  of  and 
QgQ  for  the  period  August-December  1973»  obtained, re£ 
pectively,  by  the  methods  described  in  the  1974  re- 
port and  in  the  present  one,  /s  it  car.  be  seen,  the 
variation  of  both  curves  ie  very  siniiTar  In  this  pe- 
riod, in  which  the  mentioned  distortion  wbb  not  pre- 
sent, 

3.3  I)lS0Ubi3lo]:. 

In  fig.  12  besides  the  running  mean  curve,  we 

shov'  the  values  of  the  exospheric  temperature,  T®  , 
used  in  this  report,  with  its  riuining  mean  of  27  days 
and  the  mean  doily  solar  flux  in  10,7  cm,  ^ as 
well  as  its  rimning  mean  for  three  soler  rotations, 

^10,7  • 

As  con  be  scon,  besides  the  variation  due  to  the  so- 
ler rotation,  the  daily  solar  flux  has  another  small 
varletlon  during  this  period,  v/ith  the  winter  values 
slightly  lower  than  those  of  numuer.  This  last  small 
variation  is  better  appreciated  in  the  running  mean 
curve. 

The  variation  of  the  exospheric  tenpereture  and  its 
ruiuilng  mean  is,  as  expected,  similar  to  the  variation 
of  the  solar  flux. 

About  the  Qgo  curve,  perhaps  the  difference  between 
the  two  maxima  could  be  related  to  the  difference  of 
solar  activity  of  the  two  periods.  To  the  same  cause 
could  be  attributed  the  difference  between  the  va- 
lues of  August  end  January,  although  such  differen- 


Running  mt-an  of  the  C90  and  Cq  values  obtained  by  two  different  procedures 
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cers  ftj’o  f'C  Mi.' 1 ivt  co Faloti  oen  be  drawn 

from  them.  But  v/hat  aecms  clear  in  that  neither  the 
solar  flux  nor  the  exospheric  temper*  ture  ohow  a 
aeoBonal  variation  tliat  can  explein  the  kind  of  se- 
miannual variation  of  the  Q^q  values. 

It  seene  tliat  we  are  left  with  the  onl;/  explanation 
that  the  shape  of  the  iigg  curve  should  be  attributed 
to  chanf^es  in  the  atmospheric  coir.pnr.j tion,  with  an 
increase  of  the  fO]  / [No]  concentration  relation  in 
the  mfixima  of  the  curve, 

2he  experiDiental  ciata  of  otmonydierif'  composition  seems 
to  agree  with  ouch  variations  v/ith  maxima  in  equi- 
noxes and  minima  in  suif'mer  and  winter  (liSyr  and  Maha- 
Jan  (1971),  Von  r.ahn  (1972),  Jacchia  (197^)),  Aa  aen- 
tiouod  in  a previous  report  (Alberc*  and  Gald6n,1974) 
the  annual  variation  can  be  caused  by  transport  of 
atomic  oxigen  from  the  surauer  to  the  winter  hemisphe- 
re caused  by  the  heating  of  the  atmosphere  of  the  sun 
mer  hemicphere.  The  semiannual  variation  can  be  pro- 
duced by  a similar  mecheniom,  v’lth  a hc*'ting  source 
in  the  auroral  zones,  that  produces  a displacement  of 
air  rich  in  atomic  oxigen  towards  lower  latitudes 
(Mayr  and  Volland  1971  and  1972),  This  heating  sorce 
can  be  related  to  the  semionnuol  variation  of  geomag- 
netic activity.  In  fact,  Strickland  and  Thoma8(1976) 
have  foiuid  evidence  of  transport  of  0 from  high  to 
lov;  latitudes,  at  hlght  altitudes,  during  magnetic 
atormo,  and  Prhls  and  Pricke  (197<^)  have  deduced, from 
data  of  ESHU  4,  enhancement  of  the  N2  denalty  in  re- 
gions of  high  and  middle  latitudes  duriiv:  periods  of 
increasing  magnetic  activity. 

An  increase  of  the  [0]  / fN2]  relation,  shell  alto 
produce  an  increase  in  the  density  of  Icnization.Now, 
in  the  Observatorlo  del  Pbro,  the  maxima  of  the  to- 
tal  electron  content,  always  have  been  foiuid  at  aqu^ 
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noxeo  (GaldAn  1968,  Gpld6n  and  Alberoa(1970))  and, 
with  the  exception  of  a year  of  very  high  aolar  a£ 
tivity,  the  aumiaer  values  have  always  been  higher 
than  the  winter  onea,  iiuchroBUlts  seema  to  support, 
the  idea  tliat  the  einall  winter  anocialy  of  the  va- 
luen  of  Q90  nay  not  be  eignif leant.  We  may  also  no- 
tice that,  from  the  stations  analyzed  by  Kosttr, on- 
ly Auckland  has  a semiannual  variation  of  Q bigger 
than  the  annual  one  (as  is  the  case  of  Ebro),  and 
it  is  also  the  station  that  do  not  present  vflnter 
anomaly  in  the  values  of  the  total  electron  content, 
another  similarity  with  Ebro, 

The  results  of  Titheridge  (1974)  suggest  a latitudj^ 
nal  variation  of  the  amplitud  of  the  annual  and  se- 
miannual variation  of  Q , The  results  of  Ebro  seem 
to  confirm  this  suggestions  as  far  as  the  semian- 
nual variation  is  concerned,  but  do  not  agree  for 
the  annual  one.  Nevertheless,  it  is  to  be  noticed 
the  great  difference  of  longitude  between  the  sta- 
tions analyzed  by  Titheridge  and  the  Observatorio 
del  Ebro,  If  there  la  a longitudinal  variation  of 
QgQ,  it  la  still  possible  that  a latitudinal  depen- 
dence of  the  amplitude  of  the  annual  variation  oould 
be  compatible  with  the  Ebro  results, 

4,  CONCLUSIONS, 

A method  hae  been  developed  to  obtain  daily  values 
of  the  Integrated  ionization  rate  of  the  atomic  oxl- 
gen,  Q90t  from  total  electron  content  data. 

The  method  has  been  applied  to  the  data  recorded  at 
the  Obaervatcrlo  del  i^ro  dxiring  the  period  August 
1973-July  1974  and  dally  values  of  Q90  for  this  pe- 
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rlod  have  been  obtelned. 

Q — ? —1 

The  mean  value  of  them  la  *3  1,0  10  cm  s and 
the  daily  values  show  a semiannual  variation  v/ith 
maxima  on  October  and  on  February-liarch  and  minima 
on  January  (the  sharpest  one)  and  on  summer.  The 
October  maximum  is  slightly  higger  than  the  Pebrua 
ry-March  one  and  the  minimiun  of  summer  io  slightly 
lovver  than  the  winter  one. 

The  difference  between  the  two  maxima  could  be  due 
to  a variation  of  the  solar  activity  that,  during 
the  analyzed  period  shows  a minimum  in  February  and 
klarch.  To  the  same  cause  could  be  attributed  the  dif 
ference  between  the  Qtjo  values  of  August  and  Janua- 
ry, However,  the  differences  in  both,  the  solar  ac- 
tivity and  the  Q90  values,  are  too  small  to  draw  fi- 
nal conclusions. 

An  harmonic  analysis  of  the  O90  values  gives  a first 
harmonic  of  15'»  amplitude  of  the  mean  value,  with 
maximum  on  December  the  1st  and  ninimtnn  on  June  the 
1st,  The  amplitude  of  the  second  harmonic  is  20^ 
with  maxima  on  September  the  28th  and  March  the 
29th-30th,  and  minima  on  December  the  2Bth  and  June 
29th, 

These  temporal  variations  of  Qgo*  seems  to  be  caused 
by  changes  in  atmospheric  composition,  since  cannot 
be  explained  by  parallel  variation  of  other  geophysl 
cal  parameters.  These  variations  indicate  an  increa- 
se of  the  concentration  relation  [0]  / [N2]  in  the 
maxima  of  the  Q90  curve. 

The  comparison  with  results  of  other  authors  agrees 
with  a latitudinal  variation  of  the  meon  value  of 
QgOt  with  maximum  towards  the  equator  and  also  with 
the  Q90  variation  with  the  solar  activity. 

Data  from  more  stations  at  different  latitudes  are 


needed  to  establish  a latitudinal  dependence  (if  any) 
of  the  Q90  se'^Bonal  variations,  i/ith  the  few  data  we 
have,  it  nay  be  said  tentatevily  that,  al  middle  la- 
titudes, the  amplitude  of  the  aemiannual  variation  in 
creases  with  latitude.  On  the  other  hand,  the  results 
of  Ebro  taken  in  conjunction  with  the  results  of  other 
stations,  do  not  seem  to  support  the  idea  of  a regular 
latitudinal  dependence  of  the  annual  variation,  unless 
there  is  also  a longibudinal  effect. 
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append  Pi. 

PARAi'IETEN  . OP  TH-:  ATIIOOPJIL'NI'.!  MODEL, 

The  lower  limit  has  been  fixed  to  a height  of  Zq  = 
= 120  Km,  and  the  data  at  this  level  are  indicated 
by  the  sub index  0, 

If  Tob  is  the  exospheric  temperature,  the  tempera- 
ture profile  with  height  is  given  by 

T (A-1) 


and  the  niunerical  density  of  the  j component  by 
(ofr,  2,2  for  its  deduction 


n.  =n. 

I /• 


T. 

T 


» / * 
^ -u— 


H 


where 


H.j  =IKT„)/(mjC)  and  y.=U1/lbH.j ) 

with  K > Bol?iman  constant  and  mj  molecular  weight 
of  the  j element. 

From  the  six  parameters  of  (A-1)  and  (A-2)fOne  is 
• oonstantf  another  two  depend  only  on  T«  and  the 

other  tkroe  are  function  of  T«  and  mj.  The  value 
of  T«  la  obtained  by  the  method  given  by  Jaoohla 
(1971). 
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The  formulae  for  its  obtentlon  are  : 

C:910 

b = 0,028 -T.  10"  ® 

T,s212.4  ♦1.(0,20994 -0.6  10"^  (T.  -1» 

s.  =( 8,42 ♦ T.( 2.75  10" ® (T.  -l)-7.52225  10‘^  )♦ 

♦(m./16-l)|4.822  ♦ 1.(2,298 5 10"®|T,  -1)-5.0127  10‘^). 

-0.5  10'®(T.  -1)^)10'^ 

k.= 

i 

for  the  0 : (1,25  ♦1.5  10’^  T.  ) 10’^ 

" Nj  : (0.83^T.(2.505  10‘’^^0.5  10'®  (T. -1)))10‘^ 

••  ••  02  ; (-8.76^T.(3.41375  !0'2  -3.9154  10“®  (T. -1  ))♦ 

♦11/3)  4.6  D*®  (T.  -1)®)  10‘^ 

log  n.  = 

for  the  0 : II.1212-T.(I.7508810’^-2.4465  10'^(T.-1))- 

-(7/6)  10'’°  (T.  -D® 

*'  '•  N2  : 11.0755^T.(1.16895  10"® -1.2137  10"® (T. -!))♦ 

♦(13/3)10"’°  (T. -D® 

'*  ••  O2  : 10.1168  ♦7.(1,45949  10"® -1,51335  10"® (T. -!))♦ 

♦ 5.5  10"’°(T.-1)® 

The  temperatures  are  given  in  degrees  Kelvinf  the 
heights  in  Km  and  the  densities  In  particles  per  om?. 
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U = 32  TdllP)  = 700,0  *K 
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= 20  T(IIIP)  = 700. 0*K 
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1.9999 

l.<"<-’09 

- 2 

600 

700.0 

1.6562 

1.6481 

- 2 

620 

700.0 

1.3146 

1.5078 

- 2 

640 

700.0 

0.9750 

0,9700 

- 1 

660 

700.0 

0.6373 

0.5345 

- 1 

680 

700.0 

0.3016 

0. 3011 

- 0 

- 49 


LI  ^ 16  T(Il:i-')  = 700.C)"'< 


z 

i' 

Lt'G  H 

LuO  ]: 

150 

50'3.0 

1U,3'’"0 

10, 35 3A 

1 

160 

540,3 

10.1721 

10.1773 

1 

180 

595.1 

9.8692 

9.8698 

0 

200 

631.1 

9.6007 

9.5986 

- 1 

220 

654.7 

9.3505 

9.3487 

- 1 

2'10 

6'/0,2 

9.1117 

9„1118 

0 

260 

680.4 

8.8306 

8.8i332 

1 

200 

687.2 

8.6551 

0.6603 

1 

3c.O 

691.6 

8.4339 

O.AAIO 

2 

.320 

694.5 

8.2159 

S.2245 

9 

•"-10 

696.4 

8.0005 

8.0100 

2 

JuO 

697.6 

7.7374 

7.  7^' 72 

2 

300 

698.4 

7.5761 

7.5856 

0 

400 

699.0 

7.3667 

7.3753 

0 

420 

699.3 

7.1588 

7.166? 

9 

440 

699. 6 

6.9523 

6.9582 

i 

460 

699.7 

6.7473 

6.7514 

1 

400 

699.8 

6.5437 

6.5A5S 

1 

50(^ 

699.9 

6,3413 

6.3415 

0 

520 

689.9 

6.1402 

6.1384 

- 1 

540 

700.0 

5.9403 

5.  ^''367 

- 1 

56C 

700,0 

5.7417 

5.7363 

- 1 

500 

700,0 

5.5442 

5.5374 

n 

•"  c. 

600 

700.0 

5.3479 

5 . 33^  '9 

- 2 

620 

700.0 

5.1528 

5.1438 

n 

640 

700,0 

4.9588 

4.9493 

- 2 

660 

700. 0 

4.7659 

4.7562 

_ 2 

680 

700,0 

4.5742 

4,5646 

- 2 

700 

700.0 

4.3336 

4.3744 

- 2 

720 

700.0 

4.1940 

4.1857 

_ 9 

740 

700.0 

4,0055 

3.9983 

_ 9 

760 

700.0 

5.8181 

3.8122  • 

- i 

780 

700.0 

3,6318 

3.6274 

- 1 

000 

700.  0 

3.4464 

5.  A 4-",  Cl 

- 1 

020 

700.0 

3.2622 

3.2614 

- 0 

840 

700.0 

3.0790 

3.0800 

0 

360 

700.0 

2.8967 

2.  8(^06 

1 

880 

700,0 

2.7155 

2.7201 

1 

900 

700.0 

2.5553 

2.5A15 

0 

c. 

920 

700,0 

2.3561 

2. 3635 

2 

940 

700.0 

2.1779 

2.1865 

0 

£_ 

960 

700.0 

2.0006 

2.0100 

2 

980 

700,0  ■ 

1.8243 

1.8340 

0 

C. 

1000 

700.0 

1.6490 

1.6506 

1050 

700.0 

1.2149 

1.2220 

2 

1100 

700.0 

0.7865 

0. 7875 

0 

1150 

700.0 

0.3639 

0.3546 

_ 9 

50 


u = 32  T(ii:?)  = 300. 


z 

T 

Log  i:  Jj.c 

L(;G  I: 

150 

548.7 

9.3798 

9.3743 

1 

160 

594,2 

9.0775 

9 . 074  5 

— 

1 

180 

662;  1 

8.5527 

8.5464 

2 

200 

707.5 

8.0844 

8.0765 

2 

220 

738.0 

7,6471 

7,6407 

2 

240 

758.5 

7.2294 

7.2262 

— 

i 

260 

772.2 

6.8252 

6.8256 

0 

280 

781.3 

6.4307 

6.4542 

1 

300 

787.5 

6.0435 

6.0493 

1 

320 

791.6 

5.6620 

5.6693 

o 

340 

794.4 

5.2850 

5.2030 

2 

360 

796.2 

4.9119 

4.9197 

2 

380 

797.5 

4.5422 

4.5492 

2 

400 

798.3 

4.1755 

4.1813 

1 

420 

798.9 

3.0116 

3.8158 

1 

440 

799.2 

3.4503 

3.4527 

1 

460 

799.5 

3.0915 

3.0923 

0 

480 

799.7 

2.7350 

2,7344 

0 

500 

799.8 

2,3808 

2.3792 

•• 

0 

520 

799.9 

2.0289 

2.0267 

1 

540 

799.9 

1.6791 

1.6770 

1 

560 

799.9 

1.3;'14 

1.3300 

0 

580 

800.0 

0.9858 

0.9856 

0 

600 

000.0 

0.6423 

0. 6439 

0 

620 

800.0 

0.3008 

0.5047 

1 

A — 


51 


M = 28  i'ClHiO  = 000.0‘’« 


■i. 

T 

LOG  II  J.X 

LuG  if 

150 

543.7 

10.5656 

l;-'.36ul 

1 

160 

594.2 

10.0964 

10.0935 

1 

180 

662.1 

9. 6315 

9.6257 

— 

1 

20U 

707.5 

9.2185 

9.2111 

— 

2 

220 

738.0 

8.8338 

8.8277 

1 

240 

758.5 

8. 4669 

8.4637 

— 

1 

26u 

772.2 

8.1122 

8.1125 

0 

280 

781.3 

7.7662 

7.7694 

1 

300 

787.5 

7.4268 

7.4323 

1 

320 

791.6 

7.0925 

7.0995 

2 

340 

794.4 

6.7623 

6.7699 

2 

360 

796.2 

6.4356 

6.4431 

2 

380 

797.5 

6.1110 

6.1186 

2 

400 

793.3 

5.7907 

5.7962 

1 

420 

798.9 

5.4721 

5. 4758 

1 

440 

799.2 

5.1557 

5.1575 

1 

460 

799.5 

4.3415 

4.8/113 

0 

480 

799.7 

4.5294 

4.8272 

1 

500 

759.8 

4.2193 

4.2152 

— 

1 

520 

799.9 

3.9112 

3.9055 

— 

1 

540 

799.9 

5.6049 

3.5932 

— 

2 

560 

795.9 

3.3006 

3.2931 

— 

n 

iL 

580 

800.0 

2,9980 

2.9903 

— 

2 

600 

■ 800. 0 

2.6973 

2.6098 

•• 

2 

620 

300.0 

2.3983 

2.3^16 

2 

640 

800.0 

2.1011 

2.0955 

•• 

1 

660 

800.0 

1.8056 

1.8016 

1 

680 

800.0 

1.5118 

1.5096 

— 

1 

700 

800.0 

1.2197 

1.2195 

0 

720 

800.0 

0.9293 

0.9511 

1 

740 

800.0 

0.64u5 

0.6443 

1 

760 

800.0 

0.3533 

0.3589 

1 

780 

800.0 

0.0678 

0.0749 

2 

-sa- 

il = 16  T(INP)  = 800, 


r» 

iJ 

T 

LOO  II  Ji.G 

LOO  IT 

/*» 

150 

548.7 

10.3633 

10.3647 

0 

160 

594.2 

10.1939 

lO.lf'78 

1 

180 

662.1 

9.9099 

9. <^110 

0 

200 

707.5 

9.6635 

9.6624 

- 0 

220 

738.0 

9.4372 

9.4361 

- 0 

240 

758.5 

9.2234 

9.2237 

0 

260 

772.2 

9.0178 

9.0203 

1 

280 

701.3 

8.8103 

0.8228 

1 

500 

787.5 

8.6230 

8.6294 

2 

520 

791.6 

8.4311 

8.4309 

2 

540 

794.4 

8;&418 

8.2505 

2 

560 

796.2 

8.0547 

8.0630 

2 

500 

797.5 

7.8694 

7. 0784 

2 

400 

798.3 

7.6858 

7.6945 

2 

420 

790.9 

7.5056 

7.5112 

2 

440 

799.2 

7.3227 

7.3291 

2 

460 

799.5 

7.1432 

7. 1^80 

1 

480 

799.7 

6.9648 

6.9681 

1 

500 

799.8 

6.7076 

6.7B91 

0 

520 

799,9 

6, 6116 

6.6113 

- 0 

540 

799.9 

6,4366 

6.4345 

- 1 

560 

799.9 

6,2620 

6,23'.'0 

- 1 

580 

800,0 

6.0099 

6,0846 

- 1 

600 

800.0 

5. 9101 

5.9114 

- 2 

620 

800.0 

5.7474 

5.7395 

- 2 

640 

800.0 

5.5776 

5.5609 

- 2 

660 

800,0 

5.4088 

5.3996 

- 2 

680 

800,0 

5.2410 

5.2515 

- 2 

700 

800.0 

5.0741 

5.C647 

- 2 

720 

800,0 

4.9083 

4.8901 

- 2 

740 

800.0 

4.7433 

4.7340 

_ 2 

760 

800.0 

4.5793 

4.5717 

- 2 

780 

800,0 

4.4162 

4.4097 

- 2 

800 

800,0 

4.2541 

4.2409 

- 1 

820 

800,0 

4.0929 

4.0891 

- 1 

840 

800.0 

3.9325  . 

3.9304 

- 1 

860 

800.0 

3.7731  . 

3.7726 

- 0 

880 

800.0 

3.6145 

3.6158 

0 

900 

800,0 

3.4568 

5.4598 

1 

920 

800.0 

3.3000 

3. ’<046 

1 

940 

800.0 

3.1440 

3.1502 

2 

960 

800.0 

2.9889 

2.c'964 

0 

c. 

980 

800.0  ■ 

2.8347 

2.8433 

2 

1000 

800.0 

2,6013 

2.6907 

n 

c. 

1050 

800.0 

2.3014 

2.3115 

2 

1100 

800,0 

1.9266 

1.9348 

2 

1150 

800,  0 

1.5560 

1,5600 

1 

1200 

800.  0 

1.1919 

1.1866 

- 1 

1250 

800,0 

0.8310 

O.8I42 

- 4 

1300 

800,0 

0,4764 

0.4425 

- 0 

1350 

800,0 

0,1256 

0. 0711 

- 12 

I 
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U = 32  T(IIIP)  = 900. 0®K 


z 

S 

LOO  If  J lO 

LOO  If 

150 

590.6 

9.4282 

n./fisA 

2 

160 

644.1 

8.1402 

9.1 -'.53 

— 

1 

180 

725.0 

8.6478 

8.6425 

* 

1 

200 

780.3 

8. 21*^6 

2).  2()06 

— 

1 

220 

818.2 

7.8166 

7.8117 

— 

1 

240 

844.0 

7.4385 

7.A';60 

— 

1 

260 

861.7 

7.0746 

7.0748 

0 

280 

073.8 

6.7207 

6.7255 

1 

300 

882.1 

6.3742 

6.578'7 

1 

320 

087.8 

6.0334 

6.0591 

1 

340 

891.6 

5.6971 

'^.7035 

2 

360 

804.3 

5.3646 

5.5707 

0 

4. 

380 

896.1 

5.0355 

5.0A06 

1 

400 

897.3 

4.7089 

4.7150 

1 

420 

098.2 

4.5850 

4.3877 

1 

440 

098.8 

4. 0655 

4.0645 

0 

460 

899.  2 

3.7443 

3.7435 

— 

0 

480 

099.4 

3.4272 

3.424  8 

* 

1 

500 

009.6 

3.1122 

3.1083 

— 

1 

520 

899.6 

2.7092 

2.7942 

— 

1 

540 

899, p 

2.4882 

2.482A 

— 

1 

560 

899.9 

2.1790 

2.12?o 

— 

2 

580 

899.9 

1.8718 

1.8657 

1 

600 

900,0 

1.5665 

1.8609 

— 

1 

620 

900.0 

1.2627 

1.2885 

— 

1 

640 

900.0 

0.9609 

0.9579 

1 

660 

900.0 

0.6600 

1.’.  6596 

— 

0 

680 

Joo.o 

0.5629 

( .3631 

0 

700 

900. 0 

0.0659 

U.0685 

1 

54 


r 


26 

T(IKF)  = 

goo.o'K 

Z 

T 

LOG  i'l  J..C 

LOG  li 

150 

590.6 

10.4044 

10.3‘'40 

- 2 

160 

644.1 

10.1473 

10.1417 

- 1 

180 

725.0 

9.7100 

9.7043 

- 1 

200 

780.3 

9.3202 

9.3217 

- 2 

220 

818.2 

8.9767 

B.f>713 

_ 1 

240 

844.0 

8.6442 

8. 6412 

- i 

260 

861.7 

8.3246 

0.3242 

- 0 

280 

873.8 

8.0141 

S.0161 

1 

300 

332.1 

7.7103 

7.7143 

1 

320 

387.8 

7.4116 

7.4170 

1 

340 

891.6 

7.1170 

7.1232 

2 

360 

894.3 

6,8257 

6.8321 

0 

380 

896.1 

6.5373 

6.5433 

1 

400 

897.3 

6.2514 

6, 2566 

1 

420 

893.2 

5.9678 

5. ‘'717 

1 

440 

898.8 

5.6863 

5.6!^83 

1 

460 

899.2 

5.4068 

5.4076 

0 

480 

899,4 

5.1292 

5.1234 

- 0 

500 

899.6 

4.8534 

4. 8510 

- 1 

520 

899.7 

4.5794 

4.5756 

- 1 

540 

899.8 

4.3070 

4.3021 

- 1 

56o 

899.9 

4, 0364 

4 . 03 06 

- 1 

580 

899.9 

3.7674 

3.7612 

- 2 

600 

900.0 

3.5000 

3. 4^’ 37 

_ 9 

620 

900.0 

3.2342 

3.2282 

- i 

640 

900.0 

2.9699 

2, *-'046 

- 1 

660 

900,0 

2.7072 

2.102^ 

- 1 

680 

900,0 

2.4460 

2.4430 

- 1 

700 

900.0 

2.1963 

2.1848 

- 3 

720 

900,0 

1,9282 

l.<'^282 

- 0 

74  0 

900.0 

1.6714 

1.6731 

0 

760 

900, 0 

. 1,4162 

1.41‘'8 

1 

730 

900.0 

1,1624 

1.1671 

1 

800 

900,0 

0,9100 

0.9158 

1 

820 

900.0 

0.6591 

0.6657 

2 

840 

900.0 

0.4095 

0.4165 

2 

860 

900,0 

0.1613 

0.1681 

2 

I 

\ 


m 


M = 16 


- 55  - 

T(IIIP)  = 900.  O^K 


0 

iJ 

T 

LOO  IJ  J-.O 

LOO  !•: 

150 

590.6 

10.3735 

10.3723 

- 0 

160 

644.1 

10.2092 

10.2122 

1 

180 

725.0 

9.9384 

9.<M09 

1 

200 

780.3 

9.7082 

9.7091 

0 

220 

818.2 

9.4998 

9.5007 

0 

240 

844.0 

9.3046 

9.3068 

1 

260 

661.7 

‘\1188 

9.1224 

1 

280 

873.8 

8.9391 

8.Q4''13 

1 

300 

802.1 

6.7640 

0.7706 

2 

320 

687.8 

3.5923 

8.6001 

2 

340 

891.6 

8.4232 

8.4318 

2 

360 

894.3 

6.2563 

C.2653 

2 

380 

896.1 

8.0911 

8.1002 

o 

400 

897.3 

7.9275 

7.9362 

2 

420 

898.2 

7.7653 

7.7753 

2 

440 

896.8 

7. 6043 

7.6114 

o 

460 

099.2 

7.4446 

7.4504 

i 

480 

- 899.4 

7.2859 

7. 2003 

1 

500 

699.6 

7.1263 

7.1311 

1 

520 

699.7 

6.9717 

6.9730 

0 

540 

899.8 

6.8161 

6.8157 

- 0 

560 

899.9 

6.6615 

6.6595 

- 1 

580 

899.9 

6.5073 

6.5043 

- 1 

600 

900.0 

6.3551 

6.3502 

- 1 

620 

900.0 

6.2032 

6.1f’71 

- 2 

640 

900.0 

6.0523 

6.0/'51 

_ 2 

660 

900.0 

5.9022 

5. 6943 

_ o 

“ C- 

680 

900.0 

5.7530 

5. 7''- 4 5 

— ^ 

* ^ 

700 

900.0 

5.6047 

5.5959 

- 2 

720 

900.0 

5.4573 

5.4^64 

- 2 

740 

900.0 

5.3106 

5.3020 

- 2 

760 

900.0 

5.1648 

5.1566 

O 

••  4. 

780 

900.0 

5.0199 

5.0124 

o 

800 

900.0 

4.8757 

4.6691 

o 

820 

900.0 

4.7324 

4.7269 

- 1 

840 

900,0 

4.5899  • 

4.5857 

- 1 

860 

900.0 

4.4461 

4.4454 

- 1 

880 

900.0 

4.3072 

4.3059 

- 0 

900 

900.0 

4.1670 

4.1673 

0 

920 

900.0 

4.0276 

4.0295 

1 

940 

900.0 

5.8809 

3.8925 

1 

960 

900.0  • 

3.7510 

3.7561 

1 

980 

900,0 

3.6139 

3.6204 

n 

u 

1000 

900.0 

3.4776 

3.4652 

2 

1100 

900,0 

2.0067 

2.8172 

3 

1200 

900,0 

2.1536 

2.1560 

1 

1300 

900.0 

1.5176 

1.5059 

- 3 

1400 

900.0 

0. 8980 

0,8523 

- 10 

1500 

900.0 

0.2942 

0. 2018 

- 19 
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r-  56  - 

K = 32  = 1000.  U^K 


z 

T 

Loc  If  JaO 

LOG  K 

1'^ 

150 

628.5 

9.4666 

9.4542 

5 

160 

689.7 

9.1895 

9.1 ‘241 

1 

180 

785.5 

8.7215 

8.7187 

— 

1 

200 

848.9 

B.5172 

C.5145 

1 

220 

894.6 

7.9481 

7.9461 

1 

240 

926.5 

7.6012 

7 . 6006 

* 

0 

260 

948.7 

7.26yi 

7.2702 

0 

200 

964.2 

6.9474 

6.9502 

1 

300 

975.0 

6.6554 

6.6374 

1 

320 

982.6 

6.5251 

6.3500 

1 

540 

987.9 

6.0215 

6,0266 

1 

560 

991.5 

5.7212 

5.7264 

1 

580 

994.1 

5.4242 

5.4289 

1 

400 

995.9 

5.1299 

5.1557 

1 

420 

997.1 

4.8-;B0 

4.8407 

1 

440 

998.0 

4.5484 

4.549.8 

0 

46o 

998.6 

A. 2608 

4.2608 

0 

480 

999.0 

5.9755 

5.9750 

0 

500 

999.5 

•^.6^'17 

1 

520 

999.5 

5.4099 

5.4060 

1 

540 

999.7 

5.1293 

5.1251 

1 

560 

909.8 

2.8515 

2.8462 

1 

580 

999.8 

•2.5749  ' 

2.5695 

1 

600 

999.9 

2.2999 

2.2947 

1 

620 

999.9 

2.0266 

2.0219 

•• 

1 

640 

1000.0 

1.7549 

1.7511 
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